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Abstract 

We first prove that the set of domino tilings of a fixed finite figure is a distributive lattice, 
even in the case when the figure has holes. We then give a geometrical interpretation of the 
order given by this lattice, using (not necessarily local) transformations called flips. 

This study allows us to formulate an exhaustive generation algorithm and a uniform random 
sampling algorithm. 

We finally extend these results to other types of tilings (calisson tilings, tilings with bicolored 
Wang tiles). 

1 Introduction 

In the last ten years, a lot of progress has been done about the study of tilings. Most remarkably, 
W. P. Thurston using work of J. H. Conway and J. F. Lagarias introduced the notion of 
height functions, which encode domino tilings and calisson tilings of a polygon P. 

The notion of height function appears to be a very powerful tool for the study of tilings. It 
has notably been extended by different authors jllj jl7j to study tiling algorithms for other sets of 
prototiles. 

For domino tilings, height functions induce a lattice structure on the set of tilings of a fixed 
polygon (see ^B])- Some important results are obtained from this structure: A linear time tiling 
algorithm |22], rapidly mixing Markov chains for random sampling J2] [22], computation of the 
number of necessary flips (local transformations involving two dominoes) to pass from a fixed tiling 
to another fixed tiling ^Hl ) efficient exhaustive generation of tilings [Jj [Hj • 

Dominoes are of particular importance to theoretical physicists, for whom dominoes are models 
of dimers, which are diatomic molecules (such as dihydrogen), and each tiling is seen as a possible 
state of a solid or a fluid. 
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The present paper tries to generalize previous results to figures which are not polygons, i. e. 
figures with holes. This is done by the introduction of an equilibrium function on edges of cells of 
the figure. With this tool, we prove that the set of tilings of any finite figure has a distributive 
lattice structure, of which we give a geometrical interpretation of this structure. To this end, we 
also need to introduce some structural notions, following works of J. Propp ^1] and J. C. Fournier 
[2] : The critical cycles, which induce forced components and generalized flips. 

Our approach is constructive; this allows us to exhibit algorithms to compute the objects 
introduced. As a consequence, we obtain an exhaustive generation algorithm and a uniform random 
sampling algorithm. 

We finish by proving that these ideas can be directly adapted for other types of tilings: Calisson 
tilings and tilings with bicolored Wang tiles. 

2 Figures in the plane grid 
2.1 The plane grid 

Let A be the plane grid of the Euclidean plane M^. A vertex of A is a point whose coordinates are 
both integers. 

A vertex v = {xi^yi) is a neighbour of another vertex v = {x2,y2) if l^^i — + \yi — 2/2! = 1- 
Hence, each vertex v has four neighbours v + (1, 0), v — (1, 0), v + (0, 1) and v — (0, 1) which are 
canonically called the East, West, North and South neighbour of v, respectively. An edge of A is the 
closed segment of straight line between two adjacent vertices. A cell of A is a (closed) unit square 
whose corners are vertices. Two cells are 4-n'eighbours (respectively 8-neighbours) if they share an 
edge (respectively at least a vertex). 

A directed graph G = (V, E) is symmetric if (f , v') G if and only if (f ', v) (z E for all v, v' £ V. 
In this work we deal with the symmetric directed graph (denoted by A"*") obtained from the planar 
grid A by replacing each edge vv' by two arcs {v,v') and {v',v). For an arc a = {v,v') of A"^ we 
denote by [a] its associated edge in A. 

A (directed) path P in a directed graph G = {V,E) is a sequence of vertices {vq, . . . ,Vk) such 
that (f j, fj+i) is an arc of G for every i = 0, . . . ,k — 1. We denote by E{P) the multiset of all the 
arcs used by the path P and by V{P) the multiset of its vertices. We say that G is connected if 
any two vertices of V are linked by a path. 

A path P = {vq, . . . ,Vk) with = vq is called a cycle. The cycle is elementary if Vi = Vj and 
i j imply {i,j} = {0, k}. In a plane graph one has two kinds of elementary cycles: The clockwise 
cycles and the counterclockwise ones. 

Let G = (y, E) be a symmetric directed graph. A function g : E ^ Z is skew-symmetric if 
g{v, v') = —g{v', v) for every [v, v') € E. Given any function h : V ^ 7^ we define its associated 
difference function D{h) : E ^ "Z hy D{h){v,v') = h{v') — h{v), for all {v,v') € E. Conversely, if 
G is connected, given a function g : E ^ Z which satisfies g{C) = for all cycle C Gp and a 
vertex wq of y, there exists a unique function h -.V ^ Z such that h{wo) = and D{h){a) = g{a) 
for all a (z E. 

Let E' be a multiset of arcs of G. We denote by g{E') the sum of the values g{a) over all the 
arcs a € E' (each arc a is counted according to its corresponding multiplicity in the multiset). Then 
g{E') = X^jjg^;/ g{a). For a path P, instead of g(E{P)) we use the shorthand g{P). 

We assume that cells of A are colored as a checkerboard. We thus have black cells and white 
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cells, and two cells sharing an edge have different colors. Let us define the spin function sp on the 
arcs of A^. For each arc a = {v,v'), the spin of a is noted sp{a) and given by: 

• sp{a) = 1 if an ant moving from v to v' has a white cell on its left side (and a black cell on 
its right side); 

• sp{a) = — 1 otherwise. 

For each clockwise elementary cycle C, one has sp{C) = 4:Dis{C) where Dis{C), the disequi- 
librium of C, is the difference between the number of black cells and the number of white cells 
enclosed by C. The result is true for each cycle enclosing a single cell, and each elementary cycle 
can be decomposed into such sqaure cycles. 

2.2 Figures 

A figure F of A is a 4-connected, finite union of cells of A. The unique infinite 8-connected 
component of \ F is denoted by Hqo- The other ones are called the holes of F. The set of 
all edges in the boundary of F is denoted by Eh{F). The set of edges in Hoc n F is called the 
outer-boundary of F and denoted by E(,b{F). Analogously, we denote by Vi,{F) the set of all the 
vertices on the boundary of F. 

Because of the two types of connectivity for cells, we replace (until the end of the paper) each 
vertex f of F such that each edge issued from v is on the boundary of F, by two vertices vi and 
V2, each of them connected to exactly two neighbours of v (see Figure 




Figure 1: Vertex duplication according to 4-connectivity of F and 8-connectivity of \ F. 

A figure defines a symmetric directed graph Gp = (Vf, Ep) such that Vp is the set of corners of 
cells of F (once duplication is done), and Ep is the set of arcs a such that [a] is a side of a cell of F. 
From this point of view, the clockwise and counterclockwise contours of each hole are elementary 
cycles Gp. 

For each elementary clockwise cycle of F [i. e. whose arcs are in Ep) we define Disp{G) as 
the difference between the number of black cells of F and the number of white cells of F enclosed 
byC. 

2.3 Equilibrium function 

Informally, we can say that we want to work as if F had no hole. To this end, the informal idea 
is to introduce values on edges which make holes disappear. Precisely, this is done by the use of 
equilibrium functions as defined below: 

Definition 2.1 An equilibrium function (denoted by eq) is a skew- symmetric function from Ep to 
Z such that sp(C) + eq{C) = 4:Disp{C), for every clockwise cycle G of F. 
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For figures without holes, it suffices to take eq = 0. 

Notice that, through the decomposition of cycles, a function eq is an equilibrium function if and 
only if the following conditions holds. 

• eq{C) = for each elementary cycle C around a cell of F. 

• eq{C) = —sp{C) for each cycle C which follows clockwise the boundary of a hole of F. 

More generally, for every cycle C of F, eq{C) is a number which does not depend on the chosen 
equilibrium value. We prove in Section 0] that every figure has an equilibrium function which can 
be efficiently computed. 

We also need some auxiliary functions deduced from the function eq. 

Definition 2.2 The functions eqr, t and b are defined as follows. 

1. eqr{a) = eq{a) — sp{a) for all a £ Ep. 

2. t{a) = eq{a) - sp{a) + 2 for all a £ Ef\ Eb{F). 

3. b(a) = eq{a) - sp{a) - 2 for all a£ Ef\ Eh{F). 
4- t(a) = b(a) = eq{a) + sp{a) for all a £ Eb(F). 

Note that for any arc a in Ep, t(a) — b(a) is either or 4. 

3 The lattice of tilings 

In this section we associate three classes of objects to a figure: Tilings, height functions and acyclic 
orientations. Our goal is the study of tilings, and height functions and acyclic orientations are some 
powerful tools to do this study. 

3.1 Tilings and height functions 

A domino is a figure formed by two cells sharing an edge, which is called the central axis of the 
domino. A tiling T of a figure -F is a set of dominoes included in F, with pairwise disjoint interiors 
(i. e. there is no overlap), such that the union of the tiles of T equals F (i. e. there is no gap). 
Each tiling T of a figure F is completely determined by the set of its central axis. The characteristic 
function of a tiling T defined from the set of arcs of F into Z is given by xt(o) = 1 if [a] is a central 
axis of a tile of T and XT (a) = otherwise. 

Let T be a tiling of F. The height difference in T, noted gp , is the skew-symmetric function 
defined by 

Va G Ep grio-) = eqr{a) + 2sp(a)(l — 2xt(«)) 
Then gria) € {b(a),t(a)} for every a € Ep. 

Let us define gp{a) as eq{a) + sp{a) for every a E Ep. It can be seen that for each pair {T,T') 
of tilings of F and every a G Ep, gria) — gp'ia) = 4sp(a)(xT(«) — Xt'(^))- Thus, if gp = gr', then 
T = T' . Hence, the function gp is a tool to encode the tilings. Moreover, for every arc a € E{,(F) 
and every tiling T of F we necessarily have gp{a) = gp{a). Thus gxio) does not depend on T for 
a € Eij{F). Additionally, gp{a) — gp'^a) G {—4,0,4}, for every a € Ep. 
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Proposition 3.1 Let T be a tiling of a figure F. For each cycle C of F, one has gxiC) = 0. 

This proposition is a generalization of a theorem by J. H. Conway ,4^ about tihngs of polygons. 

Proof, (sketch) It suffices to prove the result for elementary cycles since the height difference of 
each cycle is the sum of the height differences of the elementary cycles which compose it. This is 
done by induction on the number of cells of A enclosed by the cycle. 

The case of a cycle following the boundary of a hole is easily treated from the definition of 
equilibrium functions. We also verify that the proposition holds for elementary cycles of length 4 
around a cell. 

We now use induction. If we are not in one of the cases treated above, then the area enclosed 
by the cycle can be cut by a path in F, which induces two new cycles, each of them enclosing less 
cells of A than the original cycle. Thus, by the induction hypothesis, the height difference of both 
induced cycles is null, from which it is easily deduced that the height difference of the original cycle 
is null. □ 

Proposition 13. II guarantees the correctness of the definition below. 

Definition 3.2 For each tiling T , the height function induced by T (denoted by Ht ) is the function 
from the set Vp of vertices of cells of F (once necessary vertex duplications have been done) to the 
set TL of integers, defined by hT{wQ) = and D{hT) = Qt- 

We now give a characterization of height functions of tilings. 

Proposition 3.3 Let wq be a fixed vertex of Hoo- We denote by Tip the class of all the functions 
/i : — > Z satisfying the following properties: 

• h{wo) = 0; 

• for each arc a of Ep, D{h){a) € {b(a),t(a)}. 

For each tiling T, the function hp belongs to Tip. Conversely, for each h £ Tip there exists a 
tiling T such that h = hp- 

Proof. The first statement follows directly from the definition of gp. Let h belong to Tip and let 
C be a cycle around a cell. Then D{h){C) — eqr{C) = sp{C) since D{h){C) = eq{C) = 0. Clearly, 
|sp(C)| = 4. Since D{h){a) G {b(a),t(a)}, there are three arcs of C such that D{h){a) — eqr{a) = 
2sp{a) and exactly one arc satisfying D{h){a) — eqr{a) = —2sp{a). Thus, the set T of all dominoes 
whose central axis [a] is such that hd{a) — eqr{a) = —2sp{a) is a tiling T of F. The equality 
h{v) = hp{v), for each vertex v of F, is obvious by induction on the distance from wq to v. □ 

The proposition above allows one to consider each tiling as a height function. 

Lemma 3.4 For every pair of height functions h and h' and for each vertex v of F, one has 
h{v) — h'[v) = 0[4]. Moreover, h[v) — h'{v) does not depend on the chosen equilibrium function. 

Proof. Obvious by induction on the length of a shortest path from wq to v. □ 

Let < be the canonical order on functions: h <h' \i and only if h{v) < h{v') for all v E Vp. 
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Proposition 3.5 Let h and h' belong to Tip- The functions inf(h,h') and sup{h,h') belong to 
Tip. 

In the vocabulary of order theory (see for example 0,0) the above proposition can be restated 
as follows: {Hf,^) is a distributive lattice. 

Proof. Let hi = inf{h,h'). We shall prove that for every arc {v,v') G Ep, hi{v') — hi{v) G 
{h{v') — h{v),h'{v') — h'{v)} (the proof for sup{h,h') is similar). 

For the sake of contradiction, let us assume that there exists an arc a = {v,v') of F such that 
hi{v) = h{v) < h'{v) and hi{v') = h'(v') < h{v'). From Lemma|H31 one has h'{v') < h{v') — 4 and 
h{v) < h'{v) — 4. Then a := h'{v') — h{v') + h{v) — h'[v) satisfies a < —8. On the other hand, a = 
h'{v') — h'{v) — {h{v') — h{v)). Since h and h' belong to T-Lf, one obtains a > h{v, v') — t{v, v') = —4 
which contradicts the hypothesis. □ 

We define the following order on Tp: T < T' li hp < ^t'- From Proposition 13.31 (r^,<) is 
isomorphic to (T^f, <)• Thus (Fi?, <) is a distributive lattice. 

3.2 Forced components 

Definition 3.6 An elementary cycle C of Gp is critical if t(C) = 0. 

Such a cycle is strongly critical if, moreover, sp{a) = 1, for each arc a € E{C) \ Eb{F). 

We say that v and v' are critically equivalent if there is exists a critical cycle C such that 
v,v' € V{C). The equivalence classes of this equivalence relation are called the forced components 
of the figure F. 

Let Gp = {Vp, Ep) be the symmetric graph whose vertices are the forced components and where 
{U,U') is an arc of Gp if there exists v £ U and v' G U' such that {v,v') is an arc of Gp. This 
graph is called the graph of forced components of F. 

Notice that each boundary cycle of a hole of F is strongly critical. The boundary cycle of 
the outer boundary is critical if and only if F contains as many black cells as white cells. A 
strongly critical cycle can be deduced from each critical cycle by replacing each interior arc a with 
sp{a) = — 1 by a sequence a', a", a'" of three arcs of positive spins. 




Figure 2: The forced components of a figure and the graph of forced components. {Ui denotes the 
component containing the contours of the two holes) 

It is easy to see that for any cycle C of Gp, t(C) and eq{C) do not depend on the chosen 
equilibrium function. Thus, the notion of critical cycle only depends on the shape of the figure. 
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Definition 3.7 Let T he a tiling of F. The graph of T is the spanning subgraph of Gp, denoted 
by Gt = {Vf,E'), where a G E' if and only if gria) = t(a). 

By definition, every arc in Ef){F) is an arc of Gt- Moreover, a G Et \ Ef,{F) if and only if 
XT {a) = and sp{a) = 1, or XT(a) = 1 and sp{a) = — 1. 

The following proposition is the reason why we are interested in critical cycles and forced 
components. 

Proposition 3.8 Let C be a cycle ofGp. 

• If C is critical, then for every tiling T the cycle C is a cycle of Gt- Conversely, if C is a 
cycle of Gt for some tiling T, then G is critical. 

• If G is strongly critical, then for every tiling T the cycle G is a cycle of Gt which does not 
cut any tile ofT. 

Proof. If G satisfies t(C) = 0, then for every tiling T oi F one has t(C) = gT{G). Therefore 
t(a) = gT{o) for every a E E{G), whence C is a cycle of Gt- Conversely, if C is a cycle of Gt, 
then gT{a) = t(a) for every a E E{G). Thus, from Proposition 13.11 t(C) = gT{G) = and C is a 
critical cycle. 

For the second part, if G is strongly critical then t(a) > gT{o) implies sp{a) = 1, from which 
one knows that xt{o) = 0; this means that G does not cut any tile of T. Conversely, let [a] be 
an interior edge which does not cut any tile of T. By definition, gT{o) = sp{a) + eq{a). Moreover, 
since, a belongs to Ft, one has gT{o) = ■t(a) = eq{a) — sp{a) + 2 and finally ■sp{a) = 1. □ 

Corollary 3.9 If F has a strongly critical cycle {vq,vi , . . . , Vp) such that for each integer < i < p, 
an interior edge, then there exists no tiling of F. 

Proof. Let vj = {xj,yj) be the vertex of this cycle with xj + yj maximal, and, moreover, Xj 
minimal with respect to the previous condition. One necessarily has Vj-i = Vj + (—1,0), Vj-2 = 
Vj-i + (0, —1) and (moreover) Vj^i = Vj + (0, —1). Now, follow the cycle until a vertex Vj^2k such 
that fj+2fc 7^ Vj + {k, —k) (see Figure 




Figure 3: Proof of Corollarv 13.91 

Let T be a tiling of F; at least one tile of T must be cut by an edge of the path formed from 
the part of the cycle from Vj to Vj+2k- But this is impossible, from Proposition 13.81 Thus there 
exists no tiling. □ 
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From Corollary 13.91 if F can be tiled then there are three kinds of forced components: The 
component Uoo, the single components which are reduced to a single vertex, and the hole components 
which contain the contour of at least one hole. 

The following lemma establish an useful relation between the height functions and the forced 
components. 

Lemma 3.10 Let v and v' he critically equivalent vertices. For all h and h' in Tip, one has 
h{v)-h'{v) = h{v')-h'{v'). 

Proof. Let C = (uq, fi, . . . , Vp) be a critical cycle passing through v and v' . One can assume without 
loss of generality that v = vq and v' = v^- One has h{v') — h{v) = YIIIZq = 
Sto^ t(t;j, Wi+i) = h'{v') — h'{v), which yields the result. □ 

Let us choose one vertex vjj in each forced component U of Gp. From Lemma I3.1UI h < h' ii 
and only if h{vjj) < h'{vu) for all U G Vp. 
We define a distance on Hp by 

A{h,h') := \h{vu)-h'{vu)\ 

Notice that the distance satisfies the following equalities: 

A{h, h') = A{h, inf{h, h')) + A{inf{h, h'),h') 
aIh, h') = A(/i, sup{h, h')) + A{suplh, h'), h') 

3.3 Acyclic orientations and flips 

In this part we prove that the lattice of height functions of a tileable figure F is isomorphic to a 
lattice of a subclass of orientations of the graph of forced components. These lattices have been 
precisely studied by J. Propp jHj. 

Definition 3.11 A directed graph G is an orientation of F if G is an orientation of Gp such that 
\E{C) n E{G)\ = —\h{C) for all cycle C of Gp. We denote by Qp the class of all the acyclic 
orientations of F. 

Let h he a height function and let T he the corresponding tiling. Let G^ = {Vp, E) he the graph 
of strongly connected components of Gp, i. e. the directed graph defined hy {U, U') ^ E if and only 
if there exists an arc of Gp from a vertex of U to a vertex of U' . 

It is well known (see for example) that the graph of strongly components of any directed 
graph is acyclic. Hence, G^ is an acyclic orientation of G'i;'. 

Proposition 3.12 Let h he an element of Ti. Then Gh belongs to Qp. Conversely, for each 
G = (y, E) £ Qp there exists /i G such that G = Gh. 

Proof Let C be a cycle of Gp. Then = D{h){C) = D{h){E{G) f\ E) + D{h){E{C) \ E). 
Since D{h){E{G) r\ E) = t{E{G) n E) and, for any arc a of Ep, t(a) = b(a) + 4, one obtains: 
D{h){E{C) CiE) = h{E{C) n E) + 4|^(C) n E\. Moreover, from the definition of Gh one obtains 
D{h){E{C) \E)= h{E{C) \ E). Finally, = 4:\E{G) ^E\ + h{C). 
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Conversely, let G = {V, E) belong to Gp- Let g be the function defined by g{a) = t(a) if 
a e E and g{a) = b(a) if a G Ep \ E. We prove that g{C) = for all cycle C. Clearly 
g{C) = g{E{C) n E) + g{E{C) \ E). By definition of g, one has g{E{C) n E) = t{E{C) n E) 
and g{E{C) \ E) = h{E{C) \ E). Since for any arc a of Ep t(a) = b(a) + 4, one obtains that 
g{E{C)r\E) = h{E{C)r)E)+A\E{C)r\E\. Then g{C) = h{E{C)nE)+A\E{C)r\E\+h{E{C)\E) = 
b(C) + 'i\E{C) n £;| = 0. Thus there exists heU such that h{wQ) = and D{h) = g. □ 

Definition 3.13 Let G = iV^E) belong to Qp and let U ^ Uoo be in V without incoming (resp. 
outgoing) arcs. The graph obtained from G by an upward (resp. downward^ flip in U is the acyclic 
directed graph G^ = {V,E~^) (resp. G^ = {V,E~) ) where 

E+ =E\ {{U, U') : ([/, U') G ^} U {([/', U) : {U, U') G E} 

and 

E- =E\ {([/', U) : ([/', [/) G ^} U {([/, U') : {U', U) G E} 

An upward or downward flip in U corresponds to reversing all the arcs incident to U. If [/ is 
reduced to a single vertex v, the flip is said local. If U contains the contour of a hole, we say that 
it is a hole flip. 

Proposition 3.14 For any G = {V, E) G Qp, G^ = {V, E') belongs to Gp. 

Proof. Let C be a cycle of Gp. Since the number of arcs {U, U') used by G is equal to the number 
of arcs ([/', U) used by C, one has \E{G) r\ E\ = \E{G) r\E'\. □ 

Let us denote by h^ the unique function in Hp such that Gf^u = G^ . We have h^ = hi'n.V\U 
and \h^ — /i| = 4 in [/. The two corresponding tilings differ only around U . 

In particular, va.Tp a local flip in ?7 = {v} (see Figure |3J) is the replacement in T of the pair of 
dominoes which cover the 2x2 square centered in v by the other pair which can cover the same 
square. 




Figure 4: A local flip 

The upward flips defined above canonically induce an order on the set Qp. Given G and G' in 
Qp, we say that G <fiip+ G' if and only if there exists a sequence (C/q, . . . , Up-i) of vertices of Gp 
and a sequence (Go, Gi, . . . , Gp) of graphs of Qp such that Go = G, Gp = G' and, for each integer 
< i < p — 1, Gj+i is deduced form Gj by an upward flip. 

Proposition 3.15 Let h and h' belong to Tip. Then h <' h' if and only if Gh 1^fiip+ Gy. More- 
over, in this case, one can pass from Gh to Gh' by a sequence of A{h, h') / A flips. 

Proof. The direct part of the proposition is proved by induction on the quantity A{h,h'). The 
result is obvious if A{h, h') = (i e. h = h'). 
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Now, let us assume that A(/i, h) ^ and h < h' . We will prove that there exists a forced 
component U such that, for each vertex vjj of U, h[vu) < h'{vu) (which implies h{vu) < h' {vjj) — 4) 
and an upward flip can be done from Gh on U. 

Let Uq be a component such that h{vuo) < h'{vi/^). If an upward flip can be done from Gh on 
Uq, then we are done. Otherwise, there exists an arc {vq, vi), with vq in Uq and vi in another forced 
component Ui such that D{h){vo,vi) = h{vo,vi). From Lemma [3.101 we have: h'{vui) ~ ^iyui) = 
h'{vi) — h{vi) and h'{vuQ) — h{vuQ) = h'{vo) — h{vo). Moreover, h'{yi) — h{yi) = D{h'){vo,vi) + 
h'{vo) -h{vo) - D{h){vo,vi). Thus h'{vu^) - h{vu^) = h'{vuo) - K^Uo) + D{h'){vo,vi) - h{vo,vi), 
which yields h'{vu-^) — h{vu^) > h'{vu^^) — h{vug) > 0. 

Either an upward flip can be performed on Ui, or the same argument can be repeated from Ui 
to obtain an arc (fi,f2) from Ui to another forced component U2, D{h){vi,V2) = h{vi,V2)- By 
repeating the process, there are two possibilities: Either one obtains a forced component Ui on 
which no upward flip can be done, or an infinite sequence (JJi)i,z^ is obtained. But since F is finite, 
the second possibility would imply that there exists a finite subsequence {Uj, Uj+i, . . . , Uk) such 
that Uj = Uk, which is a contradiction since Gh is acyclic. 

We have now proved the existence of a forced component U such that an upward flip can be 
done from Gh around U to obtain a function /i^. Notice that IS.{hF , h') = A(/i, h') — 4. This proves 
that by induction one can pass from h to h' with A(/i, h')/A flips. 

The direct part of the first part of the proposition is obvious. □ 

Corollary 3.16 The function ip : Tip Qp defined by ip{h) = Gh is an order isomorphism between 
{Hf,<) and {QF,<fiip+)- 

Corollary 3.17 Let h and h' belong to Tip. The number of successive flips needed to pass from 
Gh to Gh' is A(/i, /i')/4. Moreover, the components U on which flips are done in such a sequence 
are those such that h{vu) ^ h'{vij). 

Proof. A flip changes A(/i, /i')/4 by one unit. Thus A(/i, /i')/4 is a lower bound and the bound 
is reached if each flip lets the quantity decrease. Thus the bound can be reached only if the 
components U on which flips are done are precisely those such that h^vjj) ^ h'{vu)- 

Conversely, we have seen that for h < h' , one can pass from Gh to Gh' by a sequence of 
A(/i, /i')/4 flips. For the general case we use inf{h,h'). Passing through Ginf(h,h')j one can pass 
from Gh to Gh' by a sequence of A(/i, m/(/i, /i'))/4 + A(in/(/i, /i'), /i')/4 = A(/i, /i')/4 successive 
flips. □ 

Corollary 3.18 For each pair of tilings {T,T'), it is possible to pass from T to T' by a sequence 
of local flips if and only if hp = hp/ for all the vertices on the boundary of F. Moreover, in this 
case the number of local flips needed is 'Yl,veEp \ ^t{v) — hpi{v)\. 

Proof. This is a special case of the previous corollary. □ 
3.3.1 Preeness and rigidity of tilings 

Informally, our tools allow one to see what is forced and what one has to choose in order to tile a 
figure. 
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Lemma 3.19 Let h belong to Tip- The function h is maximal if and only if for each forced 
component U there exists a path P from Uoo to U in Gh- 

The function h is minimal if and only if for each forced component U there exists a path P from 
U to Uoo in Gh- 

Proof. If there exists U' such that there is no path from Uoo to U' in G^, then by taking a longest 
P' which finishes in U' one clearly deduces the existence of a component U ^ Uoo of Gh with no 
incoming arcs. Thus an upward flip can be done in U and h is not maximal. 

Conversely, if /i 7^ /^max then it is possible to pass from h to /imax by a sequence of upward flips. 
The first component u on which a flip is done in such a sequence has no outgoing arc. 

The proof for /imin is similar. □ 

Proposition 3.20 For each component U such that U 7^ Uoo, one has h^[^{vu) 7^ /imax('f^(/)- 

Proof. Assume that hminivu) = h^axivu)- Take an sequence of upward flips from /imin to /imax- 
Since no flip is done on U, no upward flip can be done on each component on a path of Gh^^^ from 
U to Uoo- Thus this path is also a path of Gh^^^- By concatenating it with a path of Gh^^^ from 
Uoo to U, a cycle appears in Gh^^^, which is a contradiction. □ 

Corollary 3.21 A vertex v ofVp belongs to Uoo if f^nd only if h{v) = h'{v) for all h, h' G Tip- 

An arc a links two vertices of the same component if and only if D{h){a) = D{h'){a), for all 
h,h' e Tip- 
Proof. If V belongs to Uoo, then h{v) — h'{v) = h{wo) — h'{wQ) = from Lemma 13.101 Otherwise, 
hmin{v) 7^ /imax('f^) according to the previous Proposition. 

If an arc a links two vertices of the same component, then D{h){a) = D{h'){a) from Lemma 
I3.1fll Conversely, if an arc a does not link two vertices of the same component, then one of its 
vertices is in a component U such that U 7^ Uoo- Take an sequence of upwards flips from /imin to 
^max- A flip of the sequence is done on U , which implies that there exists a pair (/i, h') of functions 
such that D{h){a) ^ D{h'){a)- □ 

4 Effective construction and algorithms 

4.1 forced components 

When F can be tiled (which is the interesting case), the graph of forced components can be 
constructed in polynomial time: Given a tiling T, we have to construct the graph of strongly 
connected components of Gp, which can be done in linear time (see j2Uj or [S]). 

We know from matching theory that there exists a 0(n^/^) algorithm ^0] to obtain such a tiling 
r, where n denotes the area of F. N. Thiant j21j gives an algorithm which is linear in the area 
enclosed hy F [i. e. the sum of the area of F and the areas of the holes). 

4.2 Construction of an equilibrium function 

An equilibrium function can be exhibited using cut lines (see also JHl) as follows (see Figure EI): 

For each hole Hi of F, we (arbitrarily) fix a vertical segment L, = [pi,]?^] (which is called a cut 
line issued from Hi) of M'^ such that pi is the central point of a highest cell of Hi; there exists a 
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positive integer n.j such that p[= Pi + {0,ni), the vertex p'- is not in F, and, for each integer such 
that < n'^ < Ui, the point pi + (0, n'-) is the central point of a cell of F. Hence, the point p'- is the 
central point of a cell of another 8-connected component, Hj of \ F, with j i (and, possibly, 
j = oo). 



Hn 



H, 



step(l) = -4 



step(2) = 4 



step(3) = 8 +step(l) +step 2 = 8 



Figure 5: Computation of an equilibrium function by step values. 

We say that Hj is the (immediate) predecessor of Hi. This construction yields a directed tree 
whose vertices are 8-connected components of \ F. This tree is rooted in i^oo- We inductively 
define the step value of a hole Hi (denoted by step{i)) as follows: Let Ci = {vifl, Vi^i, . . . , be a 
cycle which (clockwise) follows the boundary of Hi. We state: 

step{i) = ^ step{j) - sp{Ci) 

j\Hj has Hj for predecessor 

We can now define the equilibrium value of an arc {v,v') by: 

• eq{v,v') = step{i) if v' is the East neighbour of v and the line segment [v,v'] crosses the cut 
line Lf, 

• eq{v,v') = —step{i) if v' is the West neighbour of v and the line segment [v,v'] crosses the 
cut line Lf, 

• eq{v, v') = if the line segment [u, v'\ crosses no cut line. 

This function satisfies the conditions of the definition: It is obvious for cycles surrounding 
single cells, and if Ci is a cycle which clockwise surrounds the hole Hi, one has eq{Ci) = step{i) — 

^j\Hj has for predecessor '^*^^'(.?) ~ —sp{Ci). 

The equilibrium function defined above has a specific property which can be used for algorithmic 
arguments: For each pair {v,v') of vertices of Vp, there exists a path P from v to v' such that for 
each arc a of P, eq{a) = 0. This yields that there exists a spanning tree Tp of Gp such that for 
each arc a of Tp (seen as a symmetric graph), eq{a) = 0. 
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For each arc a of Ep, \eq{a)\ < 4n (this upper bound can be reached by taking a cycle which 
uses a and cuts no hne, except a). A precise study shows that such an equihbrium function can be 
constructed in 0(n log n) time units, where n denotes the number of cells of F. 

4.3 Minimal tiling 

For figures without holes, W. P. Thurston [22] has exhibited an algorithm that builds the minimal 
tiling. For the general case, this algorithm can be generalized as follows. 

Initialization: For each vertex v, the algorithm uses a variable value h{v) and two fixed values 
low{y) and sup{v). 

From previous results, we know that for all a E Eob{F) and all tiling T, gxio) = 9f{o)- Then 
the value h{v) can be computed for all the vertices on the boundary of -f^oo- If a contradiction 
appears, then stop. Otherwise, set low{v) = sup{v) := h{v) for all v on the boundary of H^o- 

We first construct a spanning tree Tp on Gp rooted in wq. Then, for all vertices v not on the 
boundary of Hoo, we set sup{v) = t{Tpv) and low{v) = h{Tpv), where Tpv denotes the unique 
path from wq to v in Tp. 

We set h{v) = h{v) for all vertices v not on the boundary of Hoo- 

The algorithm also uses a set V consisting of vertices v such that there exists a neighbour v' of 
V satisfying h{v) + top{v,v') < h{v'). This set is computed during the initialization. 
Main loop: While V is not empty: 

• Pick a vertex v in V and update h{v) by adding 4 units. If h{v) > sup{v) after updating, 
then stop (there is no tiling). 

• Update V by adding the neighbours v" of v such that h{v") + t(v",v) < h{v) and remove v 
if necessary. 

Proposition 4.1 Given a figure F formed of n cells, the above algorithm stops after at most n? 
time passages through the loop. 

Moreover, the algorithm stops with V empty if and only if there exists a tiling. In this case, 
when the algorithm stops, one has h = hp^-^ . 

Proof. First remark that, for each arc {v,v') of Ep, h{v') — h{v) = t(f,w') in Z/4Z. This is true 
during the initialization: We introduce the cycle C formed by the concatenation of (the opposite 
path of) Tpv, Tpv' and {v,v')^ we have t(C) = 0[4] and this, together with the relations between 
t and b, gives the result. Moreover, this property is preserved by the loop. 

Thus each passage through the loop makes the sum X^ugVp ~ ^(^)) decrease by at least 

4 units since, on vertex v, sup{v) — low{v) < 4:\E(Tpv)\ < 4n (where \E{Tpv)\ denotes the number 
of arcs of Tpv), so the algorithm stops after at most passages through the loop. 

When V becomes empty, each arc {v, v') of Ep satisfies the hypothesis of Proposition 13.31 Thus 
there exists a tiling T such that hp = h. Moreover h < hp^.^ (this is true during the initialization, 
and this property is preserved by the loop), so hp^-^ = h. 

If the algorithm finds a vertex v such that h{v) > sup{v), then there is no tiling since, otherwise, 
for each tiling T, hp{v) > sup{v), which is a contradiction (clearly, from the definition of sup, 
hp{v) < supiy)). This finishes the proof. □ 
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If we take into account the implementation, the algorithm's cost is at most O(n^) time units, as 
follows: 

All the values are encoded in unary numeration, which permits to add constant numbers in 
constant time. We use the equilibrium function described in part 14.21 and the spanning tree Tp is 
chosen in such a way that for each arc a of Tp, eq{a) = 0. Thus, the initialization costs O(n^) time 
units, to compute h{v) for each vertex and h{v') — h{v) for each arcs. 

Each passage through the loop costs 0(1) time units since it consists in a fixed number of 
additions of 4 units and sign tests. This gives the time complexity. 

Of course, a similar algorithm can be designed to construct the maximal tiling of F. 

4.4 Exhaustive generation 

An exhaustive generation can be done, extending ideas of |H1 to figures with holes. Let {Ui ,1/2, ■ ■ ■ , Ug) 
be a fixed total order of forced components of F (except C/qo)- We define a total order <;ea; on 
tilings of F as follows: Given two tilings T and T', we have T <iex T' if there exists an integer 
1 < i < p such that hxiv) < h^'iv) for each vertex v of Ui and hj-iv) = hj-i{v) for each vertex v of 
Uj with 1 < j < i. 

The order <iex is a linear extension of <height, i- e. given two tilings T and T' such that 
T < height T', we have T <iex T'. 

Proposition 4.2 Let T be a (non-maximal) tiling of F and let Tgucc denote the successor ofT in 
the lexicographic order; let i denote the largest integer such that an upward flip is possible in Ui. 

The tiling Tsucc is the lowest tiling (for <height) such that hT^^^^iv) = +4 for each vertex 

V of Ui, and hj^iv) = /iTs„cc(^) Z^'" each vertex v of Uj with 1 < j < i. 

Proof. Let T' be a tiling such that hx^^^^iv) = /it(?;)+4 for each vertex v of U and hxiv) = hxsucciy) 
for each vertex v of Uj with 1 < j < i. By definition, one has Tgucc ^lex T' . 

Moreover, assume that Ht = ^Tsucc Thus, by Corollarv I3.16( one can pass from T to 

sup{T, Tsucc) by a sequence of upward flips, which contradicts the definition of the integer i. □ 

This proposition enables one to generate all the tilings of F as follows: 

Initialization: Construct the graph Gforc. of forced components, the tiling Tmim and o[Tmin)- 
Output the tiling Tmin- 

The algorithm uses a variable tiling T stored in memory, which for initialization is equal to 

mm • 

Main loop: Compute the successor of T as follows: 

• Find the last component Ui on which an upward flip can be done (if no upward flip is possible, 
then stop). 

• Construct the minimal tiling T' such that for each vertex v of U*c\C/j, hT^v) = hT'{v) and 
for each vertex v oi U, /it/(u) = /it(^) + 4. 

• Replace T by T', o(T) by o{T'), output the tiling and go back to the beginning of the loop. 

The second item of the main loop can be done in 0{n?) time units using an algorithm derived 
from the algorithm of construction of the minimal tiling (it suffices to change the initialization, 
fixing appropriate value of h{v) for v in [J^--^Uj). 
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Thus, once the initiahzation is done, the maximal waiting time between two consecutive tihngs 
is O(n^) time units. The memory space is O(n^) since for each vertex v, one has to store /it(^) 
(using unary numeration). 

4.5 Uniform random sampling 

Consider the following process: Given a tiling T, find at random a forced component C and a 
direction (upwards or downwards). If a flip can be done in C according to the chosen direction, 
then make this flip; otherwise, do not change T. Trivially, this Markovian random process is ergodic 
and converges to the uniform distribution. 

Moreover, the method of "coupling from the past" ^3 can be applied since the process is 
monotonic and one has a method to construct the maximal and minimal tilings. We thus have 
a randomized algorithm to sample domino tilings uniformly at random. The space required is 
polynomial. 

It has been previously proved m] [23] that this process is rapidly mixing for figures without 
holes. We conjecture that is remains true in the general case. 

About algorithms, the reader can also easily verify that, given a pair (T, T') of tilings, one can 

• compute in linear time if one can pass from T to T' by a sequence of local flips (it suffices to 
compare Ht and h^' on the boundary of F); 

• compute in polynomial time with (low degree) a shortest path of flips to go from T to T' 
(using the cyclic orientations and the distance) and the length of such a path. 

5 Extension to other types of tilings 
5.1 Calissons 

The same study can easily be done for calisson (i. e. tiles formed by two neighbouring cells of 
the triangular lattice) tilings to get similar results. In this case, local flips are induced by the two 
tilings of hexagons formed by six triangular cells. There are only two small differences, detailed 
below: 

• One has two types of connectivity for triangular cells (3-connectivity for cells which share 
an edge, and 12-connectivity for cells which share a vertex). Thus, some vertices have to be 
duplicated or triplicated (see Figure inj . 



• For the proof of corollary 13.91 one has to consider a part of the critical cycle with vertices 
V = (x, y) such that y is maximal. 




Figure 6: Example of "triplication". 
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Figure 7: A part of a critical cycle which implies that there is no tiling 
5.2 Bicolored Wang tiles 

The case of dominoes is a particular case of tilings with Wang tiles (i. e. 1 x 1 squares with colored 
edges, see ^3] for details). They give rise to a tiling if the colors on the edges of neighbour squares 
are compatible. 

An instance of the problem of tiling by bicolored Wang tiles is given by a finite figure and a 
coloration of the edges which are on its boundary. Hence, a domino tiling is a tiling by Wang tiles 
with one red edge and three blue edges, all the edges on the boundary of the figure being blue. 

5.2.1 Eulerian orientations 

The same study can be done for tilings with "balanced Wang tiles" such that each square has 
two blue edges and two red edges. In this case, one has hdTiv,v') = sp{v,v') + eq{v,v') if the 
corresponding edge is blue and hdT{v,v') = —sp(v,v') + eq{v,v') otherwise. This is recognizable 
as the height function for Eulerian orientations of the dual lattice, called the six-vertex ice model 
by physicists Jj. This is also equivalent to the height function for three-colorings of vertices of the 
square lattice, and to alternating-sign matrices ^^l- The results are similar to those obtained for 
dominoes. 

5.2.2 Examples with finite height functions 

For the case of "odd tiles" {i. e. tiles with exactly three edges of the same color (blue or red), see 
jl3j): One has to take a height function in Z/8Z such that hdT{v,v') = sp{v,v') + eq{v,v') if the 
corresponding edge is blue and hdT{v,v') = —3sp{v,v') + eq{v,v') otherwise. With our technique 
of equilibrium value, it is easily proved that the set of the tilings of a fixed figure has a structure 
of boolean lattice (or hypercube), even if the figure has holes. 

The case of "even tiles" (i. e. tiles with an even number of blue edges and an even number 
of red edges) is very similar, with valxivjv') = {sp{v,v'),0) if the corresponding edge is blue and 
valT{v,v') = (sp{v,v'), sp{v,v')) otherwise. These values are taken in Z/4Z x Z/2Z. 
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